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1.   INTRODUCTION 

An  essential  feature  of  transonic  "axial"  compressor  stages  is  the 
contraction  of  the  flow  passage  approaching  and  through  the  blading. 
This  is  necessary  since  the  stage  pressure  ratio  is  significantly  greater 
than  unity.   In  order  to  design  blading,  or  to  calculate  inlet  conditions 
during  operation  of  a  compressor,  the  radial  distributions  of  velocity 
and  flow  angle  are  required  in  a  radial  plane  immediately  ahead  of  the 
blading.  Axisymmetry  is  assumed. 

The  purpose  of  this  report  is  to  document  calculations  performed 
for  the  "TRANSX"  compressor  in  the  Turbopropulsion  Laboratories  at  the 
Naval  Postgraduate  School.   The  various  short  programs  developed  in 
the  process  can  be  adapted  easily  to  other  geometries  and  to  different 
applications.  Extensive  use  has  been  made  of  the  Hewlett-Packard 
Model  9830A  Calculator.  This  type  of  programmable  calculator  allows 
performance  and  design  calculations  to  be  programmed  very  quickly  by 
an  individual  investigator  or  engineer.  A  collection  of  these  mini- 
programs  will  become  the  tools  for  a  complete  design  analysis,  and 
it  is  therefore  advisable  to  document  the  approach  as  it  evolves.  The 
programs  themselves  are  not  listed  here.   It  is  the  author's  opinion 
that  the  programming  is  simple  and  best  carried  out  by  the  user. 
However,  the  programs  reported  are  identified  by  log  number  and  can 
be  listed  or  taped  if  requested. 

The  style  of  the  analysis  results  from  the  repeated  use  of  the 
"total  flow  function"  which  is  reintroduced  in  Appendix  A.   In  the 
analysis  of  internal  compressible  flows,  particularly  if  they  are 
adiabatic,  the  flow  function  provides  the  means  of  replacing  density 
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in  the  equations  with  a  well -under stood  function  of  the  velocity. 
With  a  technique  for  obtaining  velocity  from  the  flow  function  worked 
out,  as  described  in  Appendix  A,  the  treatment  of  adiabatic  internal  flows 
as  being  compressible  is  little  more  elaborate  than  the  treatment  as 
incompressible,  provided  that  a  programmable  calculator  is  available. 
Furthermore,  the  extreme  precision  of  such  calculators  allows  the  use 
of  the  compressible  flow  formulation  to  very  small  velocities. 

In  Section  2  of  this  report,  results  are  given  for  a  calculation 
of  inviscid,  irrotational  but  compressible  flow  streamlines  at  the 
rotor  face  of  the  TRANSX  compressor.   Program  IPP  used  in  these  cal- 
culations was  written  for  the  IBM360  computer.  When  reprogrammed 
for  the  HP983OA  it  will  be  reported  separately.   In  Section  2,  the 
computed  distributions  of  velocity  and  flow  angle  are  well  represented 
by  a  particular  curve  fit  described  in  Appendix  B.  The  curve  fitting 
procedure  provides  the  exact  value  at  the  outer  wall,  control  of  the 
slope  at  the  wall  and  an  exact  fit  through  three  further  points.  The 
technique  is  useful  when  data  is  known  accurately  at  only  a  few 
points.  In  Section  3  the  effect  of  a  boundary  layer  on  the  outer 
wall  is  represented  by  a  blockage  factor  and  related  to  the  boundary 
layer  displacement  thickness.   In  Appendix  C  the  blockage  factor  for 
the  TRANSX  compressor  is  evaluated  from  preliminary  test  results. 

The  analytical  representation  of  the  compressible  flow  in  radial 
planes  obtained  in  Sections  2  and  3  allows  conditions  radially  along 
and  relative  to  tbs  rotor  blading  to  be  related  to  the  total  flow 
rate  and  compressor  speed.   In  Section  k   the  flow  rate  is  related  to 
the  flow  angle  at  the  rotor  tip  for  a  given  speed.  The  results 
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graphed  by  the  HP983OA  in  Figure  k   are  used  to  show  where  the  compressor 
is  operating  in  relation  to  design  incidence  conditions.   In  Section 
5,  the  distribution  of  the  relative  Mach  number  is  calculated  for  a 
given  flow  angle  at  the  tip  and  a  given  speed.   The  results  graphed 
in  Figure  5  are  needed  to  guide  the  interpretation  of  flow  field 
measurements  as  the  compressor  speed  is  increased. 

2.   JWISCID,  LRROTATIONAL  COMPBESSIBLE  FLOW  STREAMLINES 

Compressible  flow  profiles  were  calculated  using  an  IBM  36O  computer 
program  (IPP)  which  assumes  that  stream  surfaces  are  axisymmetric  and 
locally  parabolic  in  a  radial  plane.  The  position  and  slope  of  the 
streamlines  given  by  Program  IPP  agreed  closely  with  results  of  a  flow 
net  construction  for  incompressible  flow.   [The  small  variation  from 
hub  to  tip  in  the  velocity  (0.25%)  and  therefore  also  in  the  density 
accounts  for  the  agreement] .  The  calculation  was  done  for  the  design 
flow  rate  corresponding  to  an  axial  Mach  number  at  the  tip  of  0.61+6. 
It  was  therefore  assumed  that  the  pattern  of  streamlines  was  inde- 
pendent of  the  flow  rate. 

Referring  to  Figure  1,  continuity  at  station  1  gives 

w  =  J    pVcosX  2Trrdr  (l) 

r . 
which  can  be  written  as     -1 

r  =  w/(ptVtTTr0  )  = 
1 
where  subscript  t  denotes  "stagnation"  values. 


?  r  ° 

=  w/(ptVtnr/)   =J      [pV/ptVt].cos\2.(r/r0).d   (r/r.  )  (2) 


Using  the  definition  of  the  "total  flow  function",  $,  given  in 
Appendix  A,   and  writing  R  =  r/r0  ,  Eq.    (2)  becomes 


w    RJ    §cos\d(R   ) 


(2a) 


which  can  be  written  as 


w     =     Je  cosXJ    ($/fo  )(cosx/cosX<>)  d(R   ) 
R. 


(3) 


In  the  TRANSX  compressor  \n  =  0,  so  that  from  Eq.  (3), 


where 


h     =      W/? 


?  =  I  (§/*J  cosxd(R  ) 
R 


w 


(ha) 


and  ?;  is  a  function  only  of  the  inlet  geometry. 

Using  the  data  from  Program  IPP,  the  distribution  of  the  total 
flow  function  could  be  represented  analytical  by 


$  \  1 


$ 


(i  -  r)  -  al(i  - «-)      +  ^(i  -  h 


(5) 


where  the  values 


al 

- 

1. 

7986U 

nl 

= 

0. 

38582 

*1 

= 

22 

.30U8 

(5a) 


j 


were  found  (Appendix  B)  to  fit  the  data  very  closely.  This  is 
seen  in  Figure  2.   Similarly  the  flow  inclination  was  represented 
analytically  by 


\      =  a2(l  -  R2)    +  m2(l  -  R2) 


(6) 


and  the  values 

a2  =  O.9892665 

n2  =  3.8167895  (6a) 

m   =  O.U71U926 

were  found  to  fit  the  data  very  closely,  and  this  is  shown  in  Figure 
3.  The  curve  fitting  method  is  described  fully  in  Appendix  B. 

Using  Eq.  (5)  and  Eq.  (6),  the  value  of  5  can  be  obtained  by 
integration  of  Eq.  (4a).   Since  Eq.  (5)  cannot  be  solved  explicitly 
for  $/$0as  a  function  of  R,  5;  was  evaluated  for  the  TRANSX  compressor 
by  graphically  integrating  the  data  from  Program  IPP.  A  value  of  p  = 
0.708355  was  obtained. 


3.  BOUNDARY  LAYER  BLOCKAGE 


Preliminary  measurements  of  impact  pressure  ahead  of  the  rotor  have 
shown  that  a  "boundary  layer"  exists  on  the  outer  wall  extending  appro- 
ximately 0.25"  into  the  flow.  At  smaller  radii  the  stagnation  pressure 
is  nearly  constant.   In  order  to  correct  the  inviscid  streamline  calcula- 
tion for  this  deficit  a  blockage  factor  can  be  introduced  which  can  be 
evaluated  from  measurements,  or  assumed  a  priori.   If  6*  is  the  displace- 


ment thickness  and  h*/?0   «  1, 

w  =  w  +  D„Va2TTr  6 *  where  w 

m   H*  •   •  m 

is  the  measured  flow  rate  and 
w  is  the  flow  rate  for  inviscid 
flow.  Then,  using  the  definition 
in  Eq.  (2), 


W  =  W  +  §  2- 
m    *°  r 


;* 


OUTER  WILL 

\7 


i 


Writing  k^  for  the  area  blockage  factor  usually  defined  as 


w 


B      w 


m 


(7) 


then 


or. 


W  =  kg  W  +  $ 


26^ 


1    °  ro 


f, 


,     =1  -_s-  .S£ 


and  using  Eq.  (k) , 


v    =  1  -  2-5-i- 


(8) 


Hence  to  account  for  the  viscous  blockage,  k^  must  be  assigned  a 
value  or  calculated  from  data  using  Eq.  (8),  and  then,  using  Eq.  (k) 
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and  Eq.  (7), 

is  the  value  of  the  flow  function  at  the  tip  in  terms  of  the  measured 
flow  rate.   (The  boundary  layer  on  the  hub  is  very  thin  and  is  ignored 
here ) . 


h.      CALCULATION  OF  THE  FLOW  RATE  FOR  CONDITIONS  GIVEN  AT  THE  TIP 

The  above  Eq.  (9)  relates  the  flow  function  at  the  tip  to  the 
measured  flow  rate.   It  is  desirable  to  relate  the  measured  flow  rate 
at  a  given  rotor  speed  to  the  flow  angle  or  incidence  at  the  tip.  The 
"design  point"  is  taken  as  the  condition  at  which  the  flow  is  at  the 
design  incidence  angle  at  the  tip  (P-,n  =  P-,n  )  and  the  "design"  flow 
rate  is  therefore  a  function  of  the  rotor  speed. 

If  p   is  given,  using  the  notation  of  Fig.  1, 


Uo  -^-'o 


(10) 


and 


U 


V 


tanp 


10 


(11) 


T,/T  „,  where  "ref"  denotes  a  reference  temperature  (taken  as 
t'    ref  *  v 


If  9 

Tref  "  517.8'R),  then 


V 
X   =  ; =2  = 


1 


t^o 


V 


ref  260  tanp 


10 


Ly^  J 


(12) 


where  Eq.  (9)  and  Eq.  (10)  have  been  used, 


V 


ref 


^ 


2CpT 


ref 


(13) 


and 


=v 


2CpT, 


(1^) 


as  described  in  Appendix  A. 

The  total  flow  function  at  the  tip  is  given  by 
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$        =     X  (1  -  X  ) 

*o  ov  oy 


1/y-l 


(15) 


and  using  Eq.    (9) > 


W 


m 


^1-^*0 


From  the   definition  of  W     in  Eq.    (2), 


w 


=   p.  V.  ^o     W 


(16) 


Using  Eq.  (l6),  Eq.  (15)  and  introducing  reference  pressure  and 
temperature , 


^V~J  =  prefVreKmr/  ^1*  ?  T 


X  (1-X 
o    o 


i/y-i 


(17) 


In  this  equation  pref  =  pt(e/6)>  where  6  =  (pt/pref),  Pref  is  a 
reference  pressure  (taken  as  p  _  =  29.92"  Hg)  and  a  perfect  gas  is 
assumed. 

Since  Eq.  (12)  relates  the  outer  velocity  X  to  the  referred  speed, 
N/  /~q"  and  the  flow  angle  P-,^,  Eq.  (17)  gives  the  referred  flow  rate  in 
terms  of  these  two  variables  for  a  given  blockage  factor.  The  HP983OA 
was  programmed  to  graph  the  relationship  given  in  Eq.  (17).  The  result, 
for  the  blockage  factor  (0.99)  measured  in  preliminary  tests  of  the 
TRANSX  compressor  is  given  in  Fig.  k.      (The  program  is  logged  as  HP- 

Sf 001-15). 

It  should  be  observed  that  the  design  flow  angle  at  the  tip  for  the 
TRANSX  compressor  is  £  =  650.   The  design  speed  of  the  compressor  is 


30,1+60  rpm.  In  Fig.  k   choking  of  the  inlet  flow  is  seen  to  occur  at 
about  p  =  57*  at  this  speed. 

5.   CALCULATION  OF  RADIAL  PROFILES  FOR  A  GIVEN  FLOW  RATE 

For  a  given  measured  flow  rate  and  stagnation  conditions,  W  can  be 
evaluated  from  the  definition  in  Eq.  (2).  $   is  then  given  by  Eq.  (9). 

At  any  value  of  R,  $/$   is  given  by  Eq.  (5)  using  the  method  of 
iteration  given  in  Appendix  C,  so  that  $  is  known.  The  non-dimensional  1 
velocity  X  is  obtained  by  solving  the  flow  function  as  described  in 
Appendix  A.  Then  the  non-dimensional  local  peripheral  speed  is 


x    .  »    .  jU  1  \  aL  (18) 


and  relative  speed  is 


-  g-  --VFT- 


\  =  r  =  Vx  +  Xu  (vi 


The  local  relative  Mach  number  is  then  given  by 


W        W      1  /   x 

M„  =  7 =  5-   ,  =  (20) 


w  V^v"    Vt  V(^)  (1-x2) 


and  the  flow  angle  by 

P  =  tan"1  (^-)  (21) 

This  procedure  can  be  used  to  obtain  profiles  of  relative  Mach 
numbers  and  flow  angle  for  measured  test  data. 
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Alternatively,  by  specifying  the  flow  angle  at  the  tip  $   can  be 
obtained  from  Eq.  (12)  and  Eq.  (15)  and  the  radial  profiles  calculated 
in  the  same  way. 

The  latter  method  was  programmed  for  the  HP983OA  Calculator  to 
obtain  distributions  of  relative  Mach  number  at  "design"  point  of  the 
TRANSX  compressor  as  a  function  of  the  referred  speed.   (The  program 
is  logged  as  HP-Sf001-l6) .  The  results  are  shown  in  Fig.  5.   It  can 
be  seen  in  Fig.  5  that  at  the  design  speed  of  30,U60  rpm,  only  a  small 
region  of  the  flow  near  the  hub  has  a  subsonic  relative  Mach  number. 
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Appendix  A.  Velocity  from  the  Total  Flow  Function 


The  "total  flow  function"  is  the  mass  flux  divided  by  the  "limiting" 
or  "stagnation"  mass  flux.  A  discussion  of  the  flow  function  is  given 
in  Appendix  A  of  reference  2.   In  addition  to  defining  stagnation 
values  of  pressure, temperature  and  density,  a  "stagnation"  or  "limiting" 
velocity  is  introduced  from  the  energy  equation.  Thus  if 

T2 


»t  =  »  +  I 


MD 


defines  the  total  enthalpy, 

--V 


Vt   =V2ht 


A(2) 


is  the  "limiting"  velocity  obtained  in  an  adiabatic  expansion  to 
zero  temperature . 

For  a  perfect  gas,  the  equations  relating  static  and  stagnation 

properties  become  very  simple  if  the  velocity  is  made  dimensionless 

(3) 

by  dividing  by  V  .     Writing 


x    =    v/vt 


A(3) 


then 


T/T        =     1  -  X' 


p/p.      =      (1-X^) 


y/(y-D 


>/p+     =     (1-3T) 


1/(7-1) 


A(h) 
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and  the  total  flow  function,  $ ,  is  then  given  by 


$   =  pV/ptVt  =  X(l  -  3T) 


1/(7-1) 


A(5) 


In  analysing  internal  compressible  flows,  the  flow  function  in 
this  form  can  often  be  used  to  eliminate  density  from  the  equations. 
In  many  cases,  the  flow  rate  and  the  stagnation  properties  are 
known  and  Equation  A(5)  must  be  solved  to  obtain  the  velocity.   Care 
must  then  be  taken  to  recognize  the  subsonic  and  supersonic  roots  as 
shown  in  the  following  sketch  of  Eq.  A (5): 


*■  5UIERS0NIC  FLOW 


X 
M 


*m 


1.0 


1.0 
oo 


Differentiating  Eq.  A(5), 


d$/dX     =     $[l/X   -  2X/(Y   -  1)(1   -  X   )] 


A(6) 


Ik 


For  the  subsonic  root  the  procedure  is  as  follows : 

(i)  guess  X<  (y-l/y+l);  X  =  $(l+$)  is  suitable 
(ii)  calculate  $   =  X(i  -  X2) 


j.VCt-i) 


(iii)  test  for  convergence;  1$,   -  $  |  <  e?  where  $   is  the  given 
value 

(iv)   calculate  X  =  X  +  ($   -  $)/(d$/dX)  using  Eq.  A(6) 

(v)   repeat  from  (ii)  for  convergence  at  (iii) 

For  the  supersonic  root  choose  an  initial  guess  at  (i)  such  that 
X  >  (y  -  3/fy  +  l)  and  use  the  above  procedure  without  change. 
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Appendix  B.   Curve  Fitting  with  a  =  ae  +  me 


Compressor  blading  calculations  involve  the  selection  of  a  number  of 
stream  surfaces  between  which  equal  fractions  of  the  total  flow  rate 
pass.  By  referring  properties  to  conditions  at  the  outer  radius, 
distributions  of  quantities  of  order  unity  or  less  must  be  represented. 
A  function  which  has  been  found  to  be  useful  to  approximate  smooth 
distributions  calculated  only  at  a  limited  number  of  points  is  the 
following ; 


=  ae   +  me 


B(l) 


The  shape  of  a  =   CT(e)  given  by  Eq.  B(l)  depends  on  the  range  of 
e  and  particularly  on  the  value  of  n.  Differentiating, 


da/ds   =  ane11"1  +  m 


B(2) 


from  which  it  is  evident  that  there  is  a  discontinuity  in  behavior  at 
n  =  1  for  n  going  from  values  less  than  to  greater  than  unity.  This 
can  be  seen  in  the  following  sketches: 


Case  1:  rt  <  I 


Case  2:    tt.  >  I 


i.o       £ 
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In  Case  1,  the  slope  is  infinite  at  e  =  0  and  is  given  by  the  value 
of  m  as  e  becomes  large.   In  Case  2,  the  slope  is  given  by  m  at  e  =0 
and  increases  as  e  becomes  large. 

The  function  can  be  fitted  to  given  data  in  a  number  of  ways.  There 
follows  the  analysis  for  fitting  three  arbitrary  points  exactly.   (in 
application  it  is  found  that  particular  choices  of  data  points  can 
be  made  to  improve  the  overall  fit.) 

If  subscripts  1,  2  and  3  denote  three  given  data  points,  the  equations 


n 
°1  =  ae  1     we   1 


ct2  =  ae  2   +  me  B(3) 


a,  -  ae  ,   +  me 


'3 


must  be  solved  to  give  the  values  of  a,  m  and  n. 

Eliminating  m  between  pairs  of  Eq.  B(3)  and  'a'  from  the  resulting 
pair  of  equations,  n  is  given  by 


n-1 

1 


-  62  e^Ke^0^     "  e-jovj 


B(U) 


n-1      n-1 


e0     ~e^  e-iVe->a 


-i(.eoO'2   "  620*0/ 


Q  t-|\S*p         O  r^J  o 


When  n  is  known,  a  and  m  are  given  by 


e  2O-1      "     e  -1  Op 
n  ~~ n 

S2el      "  ele2 


B(5) 
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and 


m 


=     (o±  -  aein)/ei     =     (CT2   -  ae2n)/e2  B(6) 


Thus   it  is  necessary  first  to  solve  Eq.  B(U).     Writing 


L     =     [LHS   of  Eq.   A(U)]    =  [(e^^f1   -  1]/[1   -    (^/e^"'1] 


or 


L     -     (A71  -  1)/(1  -  B71)  B(7) 

where  A  =  ej/fgj  B  =  &J&2  and  7]     =     n  -  1. 

Choose  e-,  e      and  e      such  that  0<e<e<e,  then  A  <  1,  B  >  1 
and  T)   (<  0)   is  to  be  determined. 

Differentiation  of  Eq.  B(7)   gives 

dL/dn     =     dL/dTl     =      (A71   -  B71)/(l  -  B71)2  B(8) 

From  an  inspection  of  Eq.  B(7)  and  Eq.  B(8); 

(i)     If  7]  >  0  (n  >  1)    ,  L  >  0     and   (dL/dn)  <  0 

(ii)   If  T)<  0   (n  <  1)    ,   L  >  0     and   (dL/dT))  >  0 

Note  that  as  T)  -♦  »,  L  -»  0.  Hence  the  situation  for  obtaining  a 
solution  to  Eq.  B^)  is  as  shown  in  the  sketch,  where 
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So    vepCT-,      -     Snap) 

R   =   TRHS   of  Eq.    B(U)l    =  / ±-^r 

St    ve  oCT0      "     eoO-J 


B(9) 


1   Vfc3u2 


2U3' 


If  an  initial  value  of  j\  is   chosen    (   =  T)    ) }   then  an   improved  guess 
is  given  by 


T]     =     Tl- 


R-Ll 

dL 

dTl 

1 

B(10) 


where  j  =  +  1  when  f]_    = 


+ve 


ve 


The   initial  guess  should  be  made  from  an  inspection  of  the  curva- 
ture required  to  fit  the  data. 

If  we  define  the   difference   in  slopes   as   follows : 


-     <7q   ~  <t. 


C     =     "3 


ov,   -  a- 


e2    "  el 


B(ll) 


then 


C  >  0  requires     7]  >  1   (l|  >  0) 


C  <  0  requires     T]  <  1    (t]  <  0) 
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Initial  guesses  of  7]  =  -0.5  for  C  <  0  and  7]  =  2.0  for  C  >  0  have  been 
found  to  be  satisfactory 

This  procedure  was  applied  to  represent  the  distributions  of 
total  flow  function  and  flow  inclination  angle  at  the  TRANSX  compressor 
rotor  face.  Data  for  eight  streamlines  calculated  for  equal  mass 
flux  increments  were  obtained  from  IBM  360  program  IPP   .   The 
HP983OA  calculator  was  programmed  to  accept  eight  pairs  of  data 
values,  to  calculate  a,  m  and  n  by  curve  fitting  three  specified 
pairs,  to  plot  the  resulting  function  over  the  range  of  interest  and 
finally  to  plot  the  input  data.  The  results  are  shown  in  Fig.  2  and 
Fig.  3  (for  which  the  programs  are  logged  as  HP-Sf 002-1  and  HP-Sf 001-3 
respectively).   It  can  be  seen  that  the  given  data  is  very  closely 
represented  by  the  assumed  expression,  however  the  flow  function  is 
not  represented  explicitly  in  terms  of  R  and  must  be  obtained  by 
iteration.  The  procedure  is  as  follows: 

(i)  To  obtain  an  initial  guess  for  §  that  is  not  greater 
than  unity  a  power  law  is  assumed 

I  =  1  -  0.03  (1  -  Ri)2 

(ii)  Calculate 

R   =  1  -  ax  -  mx  ,  where  x  =  1  -  $  and 
$  =  $/$o 
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12    2  I 
R   -  R  I  <  e  ?  where  R  is  the 

given  value . 


?    ?  n-1 

(iv)   Calculate  the  slope  at  R  ;  dR  /dx  =  -anx    +  m 


(v)   Calculate  $  =  $  - 


El2  -  R2 

dx 


(vi)  Repeat  from  (ii)  for  convergence  at  (iii). 
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Appendix  C.  Blockage  Factor  from  Impact  Pressure  Measurements 

Radial  surveys  of  impact  pressure  at  station  1  of  the  TRANSX 
compressor  showed  that  the  impact  pressure  was  constant  at  radial 
locations  more  than  1"  from  the  outer  wall  and  the  dynamic  pressure 
was  within  k%   of  the  core  value  to  within  0.2"  from  the  wall.  Treating 
the  region  near  the  wall  as  a  boundary  layer,  the  displacement  thick- 
ness is  given  by 

y» 

*  -[  (1"f^)ay  c(1) 

where  oo  denotes  the  value  at  r  <  U.5"  and  y  is  measured  radially 
inwards  from  the  wall.  In  terms  of  the  total  flow  function,  Eq. 
C(l)  can  be  written 

yoo 

6*  -   J(l"t")dy  C(2) 

where,  using  Eq.  h{h)   and  Eq.  A(5) 

1/ 


$  =  (p/pt)yVl  -  (p/Pt)(7"l)/y  C(3) 


In  terms  of  the  difference  A  (ins.  water  col.)  =  p  -  p, 


p/pt  - 


1 
A/p  +  1 


C(U) 


Equations  C(2),  C(3)  and  C{k)   were  programmed  for  the  HP983OA  and 
the  following  table  gives  the  results  for  TRANSX  compressor  Run  8:  " 
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y  A  (ins.  water)  (l  -  $/$  ) 

^  0  4.19               .2508818 

.05  5.66                .1326629 

.1  6.1+9               .0732U42 

.2  7.27                .0211151 

.3  7.46  8.89409  x  10"3 

.5  7.57  I.898II  x  10"3 

.6  7.59  6.32237  x  10 

.7  7.58  I.26U9U  x  10"3 

.8  7.58  1.26494  x  10"3 

.9  7.60                  0 

1.0  7.60               0 

($   =  .07305442  for  y  =1.0") 
00  00 

The  integral  in  Eq.  C(2)  was  calculated  from  the  tabulated  data 

using  quadratic  integration  described  in  Appendix  B  of  reference  2. 

For  this  data  6  =  0.0210594  and  using  Eq.  8  with  5  =  .70836,  the 

blockage  factor  k^  =  1  -  2.82344  ^-  =  O.9892. 

o 
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Figure  1.   GEOMETRY  AND  PRINCIPAL  NOTATION 


2k 


0.97 
1.0  1 


0.9B    -l/l0     0.99 


0.9  ■ 
R/R0 

0.B  ■■ 
0.7  - 
0.G  - 


IRRN5DNIC  COMPRESSOR 
CRLCULRTED  INLET 
FLDWFIELD 

X   COMPUTER  PRDERRM 
CURVE  FIT 


0.S  4 


1. 00 


Figure  2 .      RRD I RL  D I 5TR I  BUT  I DN  DF 
0.0  0.2 


0.B  - 


0.E  ■■ 


0.H 

H — 


0.E  RRDIRN5 

H 


TRRN5DNIC  CDMPRE55DR 
CRLCULRTED  INLET 
FLDWFIELD 

X   COMPUTER  PROGRAM 
CURVE  FIT 


Figure  3.   RRDIRL  DISTRIBUTION  OF  THE  FLOW  RNELE 
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